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We show that invariance properties of the Lagrangian of an incommensurate system, as described 
by the Frenkel Kontorova model, imply the existence of a generalized angular momentum which is an 
integral of motion if the system remains floating. The behavior of this quantity can therefore monitor 
the character of the system as floating (when it is conserved) or locked (when it is not). We find that, 
during the dynamics, the non-linear couplings of our model cause parametric phonon excitations 
which lead to the appearance of Umklapp terms and to a sudden deviation of the generalized 
momentum from a constant value, signalling a dynamical transition from a floating to a pinned 
state. We point out that this transition is related but does not coincide with the onset of sliding 
friction which can take place when the system is still floating. 

PACS numbers:05.45.-a, 45.10.-b,45.05.+x,81.40.Pq 



I. INTRODUCTION 

Measuring friction at an atomic level has recently be- 
come experimentally possible [Q. Many studies of the 
dynamics of appropriate non-linear systems aiming at es- 
tablishing the mechanisms giving rise to energy dissipa- 
tion during the sliding of a body onto a crystalline sur- 
face have appeared in the literature |^-|^. The Frenkel- 
Kontorova (FK) model, which describes a harmonic chain 
interacting with a rigid periodic substrate, is particularly 
suitable to study the important case of an incommensu- 
rate (IC) lattice parameter ratio of the contacting sur- 
faces. The present study focuses on the effects of dis- 
commensuration on the dynamics. It should be kept in 
mind that a more realistic study of friction would require 
an extension to two dimensions. Coupling to the third 
dimension can be provided either by an ad hoc damp- 
ing term or by coupling to an elastic medium The 
ground state properties of this model have been thor- 
oughly studied 0. At a critical value Ac of the coupling 
to the external potential, the ground state of the system 
displays a structural transition (Aubry transition) from a 
floating to a pinned configuration. Below this threshold, 
the center of mass of the static system can be displaced 
on the substrate without energy costs. Therefore, one 
might expect a frictionless regime also in a dynamic sit- 
uation and a superlubric regime, where the chain would 
slide indefinitely, has been predicted for this case [|| . In a 
previous paper |^ we have pointed out that the inherent 
non-linear coupling of the center of mass (CM) motion to 
the phonons leads instead to an irreversible decay of the 
CM velocity. The essential mechanism for the transfer of 
kinetic energy from the center of mass to the internal vi- 
brations is the parametric resonant excitation of phonons 
mediated by ordinary resonances with phonons related to 
the modulating potential. 

Here we show that this type of mechanism has another 
important consequence, namely it causes the appearance 
of Umklapp terms, signalling a dynamical transition in 
the system from a floating to a pinned state. We have 



studied this phenomenon by identifying a new quantity, 
which we call generalized angular momentum (GAM), 
which is an integral of motion only if the system is in 
a floating IC phase, reflecting the invariance of the La- 
grangian of the model for a phase shift in this state. We 
show that this invariance is equivalent to the absence 
of Umklapp terms. By means of numerical simulations 
we show that the temporal behavior of the GAM is a 
powerful probe both of the (in)commensurability of the 
groundstate configuration and of the dynamical phase in 
which the system is during motion. Simulations where 
the incommensurate ground state is given an initial ve- 
locity show that the GAM remains conserved up to a 
well defined time where a sudden jump takes place. We 
have been able to relate this change of behavior from con- 
served to non conserved to the appearance of Umklapp 
terms. An important finding is that this floating-pinned 
transition does not coincide with the onset of friction. 
It was recently suggested by Popov j|] that the appear- 
ance of Umklapp terms, i.e. the conservation of quasi- 
momentum instead of momentum for crystalline systems, 
is the mechanism via which friction occurs in incommen- 
surate contacts. The present result shows that this is not 
the only mechanism. By monitoring the system via the 
GAM we can show that decay of the CM velocity may 
occur already in the floating phase. The onset of friction 
and the appearance of Umklapp terms are both caused by 
non-linear couplings and resonant phonon excitations in 
the system but remain two distinct phenomena occurring 
at different times. 

In Sec. m we describe the construction of the GAM by 
deriving it from the Lagrangian for the system in Fourier 
space a nd de fine conditions under which it is conserved. 
In Sec. Ill A we present results of numerical simulations 
which confirm the validity of our analytical derivation 
and underline the usefulness of the GAM to discriminate 
between commensurate and floating IC and pinned IC 
phases, respectively. Subsequently, we examine in Sec. 
[II B the relationship between pinning and Umklapp and 
show the presence of a well-defined transition time. In 
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Sec. [V we present conclusions and perspectives of this 
work. In the Appendix we provide the reader with an 
exphcit proof that the GAM is an integral of motion in 
absence of Umklapp. 



II. CONSTRUCTION OF A GENERALIZED 
ANGULAR MOMENTUM 

In this section we will construct a generalized angular 
momentum for the dynamical FK model, as described in 
Ref. Q . We remind the reader that this model represents 
a chain of N particles which interact with each other via 
a first-neighbor harmonic potential and are subjected to 
an external, spatially periodic, potential of strength A. 
The FK Hamiltonian reads: 
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where the u„ are the particle positions and Pn their mo- 
menta. The ratio between the modulation period of the 
external potential m and / (the equilibrium distance be- 
tween the atoms of the chain for A = 0) is taken to be 
irrational, i.e. the system is incommensurate. In our 
calculations we take m = 1 and I — t = (golden 
mean). In the numerical implementation for a finite sys- 
tem of N particles we impose periodic boundary condi- 
tions 
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It is important to notice that since wave vectors are de- 
fined modulo 27r, the Kronecker deltas in Eq. (^ should 
be read as: 
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Umklapp is present whenever this relation is satisfied 
with s 7^ 0. It is clear that the occurrence of Umk- 
lapp depends on the modes Xk which are not negligible, 
and on the choice of the (extended or reduced) Brillouin 
zone in which k is represented. It is known that in the 
groundstate, for a coupling A well below the critical value 
Ac, which for this model assumes the value Ac = 0.154..., 
the modes with wavevector nq have an amplitude which 
scales as A'"' . This number \n\ is therefore a natural label 
to represent the modes; we define n{k, q) as the smallest 
(in absolute value) number which satisfies: 



k = n{k, q)q mod (27r) 
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For a finite system with N particles, where k can be rep- 
resented in the reduced Brillouin zone as k = K{2tt/N), 
K G (-l/2iV, 1/2A^], this can be rewritten as: 
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This implies that we have to choose commensurate ap- 
proximants for the equilibrium distance I. By expressing 
I as ratio of consecutive Fibonacci numbers we obtain 
approximants which satisfy the condition I x N = M x 1 
with M and N integers. Let us introduce the modula- 
tion wavevector q = 2ttI ~ 2tt{M/N) and the position 
and momentum of the CM of the chain of atoms: 
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The equations of motion for the deviations Xn — Un — 
nl — Q from the equilibrium positions in the uncoupled 
chain are then given by: 

Xn — Xn+1 + Xn-1 - 2Xn + A COS {qn + 2'KXn + 2nQ) (4) 

As noted in Ref. , in the weak coupling regime it is use- 
ful to move to Fourier coordinates Xk = jf X^n e~*'^"a;„ 
with k = 2'kK/N . The phonon dispersion of the chain 
for A = is denoted by oJk = toik) — 2|sin(fc/2)|. The 
Lagrangian associated with Eq. (nl) in transformed space 
becomes: 
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FIG. 1. FK model for = 377, A = 0.05. (a) Phonon am- 
plitudes squared plotted as a function of the wavevector K, 
as in Eq. (^. The first two nq modes are explicitly indicated, 
(b) Same as in panel (a), relabeled according to Eq. (^. Due 
to finite numerical precision, the exponential decay with A'"' 
is apparent only for the first fifteen modes. 

In Fig. |l| we compare the phonon amplitudes for the 
ground state of the FK model for N = 377, A = 0.05, 
plotted as a function of the usual wavevector label K 
(panel (a)), as well as reordered according to the label n 
(panel (b)). Note that, due to finite numerical precision, 
the scaling behavior is hidden in numerical noise after 
the first fifteen modes. 

The use of n as a mode label makes apparent the fact 
that there is no Umklapp in the ground state of the FK 
model in the modulated phase for A < Ac. In fact Umk- 
lapp would imply the presence of a non vanishing term: 

XniqXn^q ' ' ' Xn^q ] + ^2 H = sN with S ^0 

(9) 

The joint amplitude of this term would be: 

X\ni\+-\n,„\ < X\''\N ^jQ^ 

which vanishes for N —>■ oo. 

The absence of Umklapp terms is directly related to the 
existence of a free floating phase, which is a well known 
invariance property of the FK model. In the present no- 
tation it amounts to the invariance of the Lagrangian for 
the transformation: 

Q^Q + q^/2n (11) 
Xk ^ Xke^'^'^ (12) 

This invariance is related to the existence of a zero- 
frequency Goldstone mode in the system. This mode 
is also often called phason, and should not be confused 
with the usual acoustic mode of periodic crystals. 

Having found an invariance for the Lagrangian, we can 
look for the conjugate conserved momentum. We get: 

P0 = — r = V nqX^nqXnq + T^Q = L+ —Q (13) 

The quantity represents a generalized angular momen- 
tum (GAM). It is important to realize that the invariance 
of the Lagrangian only holds in a subspace of the full 
phase space where Umklapp terms can be neglected as 
it is the case for the floating (incommensurate) ground 
state. In order to stress this point, a direct calculation of 
P0 is given in the Appendix, showing that the GAM 
is an integral of motion only if Umklapp is not present. 

This quantity is therefore an useful tool to discriminate 
between commensurate and incommensurate structures, 
and floating and locked states. In the next section we 
present numerical simulations that we carried out for var- 
ious values of the parameters of the model, showing how 
P0 is a good indicator of the phase in which the system 
under examination is. 



III. NUMERICAL RESULTS 



A. Commensurate vs. incommensurate, locked vs. 
floating 



We have performed numerical simulation in order to 
study the behavior of the GAM, as defined by Eq. ( p^ ), 
integrating by a Runge-Kutta algorithm the N Eqs. ^). 
We assign to the particles of the chain as initial conditions 
momenta Pn = Pq and positions Xn {t — 0) corresponding 
to the ground state. Fig. |^ shows simulation results for 
the same number of particles N and potential strength A, 
but for a low (r — 5/3) and a high (r — 233/144) approx- 
imant to the golden mean r, producing a commensurate 
structure and an approximate incommensurate one. The 
qualitative behavior of the momentum of the center of 
mass P is similar, whereas the behavior of Prf, in Fig. ^ 
is remarkably different being conserved only for the case 
which approximates an incommensurate system. This 
confirms that can be used as a tool to discriminate 
unambiguously between commensurate and incommen- 
surate structures. 




Time 

FIG. 2. Behavior of the CM momentum P for: (a) an in- 
commensurate configuration with A'^ = 144, r — 233/144, 
A = 0.015, Pq — 0.29; (b) a commensurate case with A'^ — 144, 
r — 5/3, A = 0.015, Po = 0.29. Note the qualitative similarity 
in the behavior of P. 

Furthermore, our numerical simulations show a re- 
markable fact. If we start the simulation with an in- 
commensurate initial condition, p^ is indeed conserved, 
but only up to a critical time tc, after which it rapidly 
deviates from its initial conserved value. This is shown in 
Fig. ^, where we can examine the behavior of and P in 



3 



a weak coupling, highly incommensurate (r = 610/377, 
A = 0.015) case. In order to check that the observed 
variation of only sets in after a critical time tc, we 
have analyzed the behavior of the quantity log{p^ — Co), 
Co being the value of at t = 0. It is evident from 
panel (c) of Fig. ^ that we can identify such a critical 
time tc where the GAM has a jump in value of various 
order of magnitude. Besides, this figure shows that, for 
t < P(ji is conserved within our numerical accuracy, 
never exceeding variation larger than 10"^*^. 



with frequency VL = 2'kP in the equations of motion of 
the system. Linear stability analysis |^,|| shows that a 
mode k grows exponentially (with rise time r) whenever 
its frequency satisfies: 

^ _ uj{k) + uj{mq - k) 
m 

for some m. The ~ symbol indicates an instability win- 
dow of relative width Wm that scales with A I™'. 
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FIG. 3. Behavior of for the parameters of the model 
as described in Fig. ^ (a) Incommensurate case: the GAM 
is constant within numerical precision, (b) Commensurate 
case: the GAM is not conserved. Note the change of scale 
going from panel (a) to panel (b). 

The critical time tc obviously depends on the coupling 
strength A < Ac and on the initial velocity Pq. We are 
currently investigating the dependence tc{\,Po), which 
turns out to be rather intricate, and plan to report on 
it in the near future. Here, we want to concentrate on 
the mechanism that causes the breakdown of the con- 
servation law, leading to a well-defined tc- We show in 
the next section that the answer lies in the appearance 
of Umklapp terms, which render the system pinned, thus 
leading to a non conservation of p^. 
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FIG. 4. TV = 377, t = 610/377, A = 0.015, Po = 0.29. (a) 
Behavior of the CM momentum P. (b) Behavior of p^. It is 
possible to see how the GAM stops being conserved, (c) In 
order to check if there is a critical time tc, we plot the quan- 
tity log{p^ — Co), where Co = P4,{t — 0). We show that tc 
can be unambiguously identified. Note that, even if, in this 
weak coupling case, the decay of the CM coincides with the 
breakdown of conservation of p^, this is not always the case. 
See also text and Fig. M. 



B. The role of Umklapp processes 

As we have pointed out in Sec. ||, is conserved only 
in the absence of Umklapp, that is to say when terms 
of the form given in Eq. have a vanishing amplitude 
for — > oo. As we have seen, this is the case for the 
groundstate of the incommensurate system. However, 
starting from the ground state, this may change during 
the dynamics due to parametric resonances. The move- 
ment of the CM with velocity P induces a modulation 



Suppose now that a mode Xk is unstable. Let n(fc, q) 
be its label according to Eq. (^. When n{k,q) — 0{N), 
this mode could lead to Umklapp as soon as its ampli- 
tude becomes finite. However, since its initial value is 
of order A'"', it requires an infinite amount of time (as 
A^ — ^ oo) to render Al"le*/'^ finite. So we have to look 
for modes k in the instability windows with n{k, q) fi- 
nite. Since the mapping n(fc, q) —> k leads to a uniform 
distribution, one indeed expects to find such a k with 
|n(fc, q)\ < 1/w, where w = Wm is the relative width 
of the joint instabilities windows. Once the amplitude 
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of the mode k starts growing (with a behavior given by 
an exponential law of the form: Ae*/^, with A bounded 
from below via the upper bound on n), also modes fc' 
with n'(k' , q) = pn(k, q) ±l,p — 2, 3, ... start to develop 
via non linear terms in the equations of motion (see Ref. 
for details), with the form: 



A(Ae*/^)f 



(15) 



Umklapp can result when repeating this process 0{N) 
times (i.e., whenpn ~ N) still gives a finite result. From 
Eq. ( [l5| ) it is clear that this will happen when Xk{t) ex- 
ceeds some threshold value, which is the case for t larger 
than the critical time tc introduced in Sec. Ill A. 
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FIG. 5. iV = 377, T = 610/377,A = 0.05, Po = 0.15. The 
appearance of Umklapp terms is shown, (a) Phonon ampli- 
tudes squared at t = 0. (b) Same at t — 1000. The mode 
K = 18 {n = 47) is starting to grow; (c) Same at t — 1500. 
The modes with n x 47, with n = 1, 2, 3 are all present, but 
there is still no Umklapp. (d) The mode with K = \9> has 
become order unity and Umklapp appears, signalled by the 
excited modes at the zone boundary, (e) Time evolution of 
log{p^ — Co). The appearance of Umklapp in panel (e) cor- 
responds to the breakdown of the conservation of the GAM. 
See also text for a detailed explanation of the Umklapp mech- 



Fig. ^ shows this mechanism at work: here we have 
taken Pq = 0.15, so that ^ = 2-kP = and N = 377 

(r = 610/377). For this value of Vl Eq. (|l|) has ap- 
proximate solutions for n = 2 for fc = K{2tt/N), with 
K — 18, 19. The corresponding values for n{k, q) are 
in the case n(18, (7) = 47 and n(19,g) — 97. Panel (a) 
shows that the phonon amplitudes at t = decay indeed 
exponentially with n. Panel (b) shows that the unstable 
mode K = \^ with the lowest value of n starts to grow 
'Ait — 1000. Modes at n = 2 x 47 are also present, due to 
non linear terms, but quadratically smaller, as explained 
above. At t — 1500 modes at n = 3 x 47 become visible, 
as shown in panel (c). There is still no Umklapp, as can 
be seen from the absence of an amplitude at the zone 
boundary. Such a term, corresponding to n = 4 x 47, 
finally appears in panel (d), aX, t = 2000. Indeed, this 
is also precisely the time at which l^igp becomes order 
unity and stops being conserved (see panel (e)). 

The mechanism described above shows that the ap- 
pearance of Umklapp terms causes a sudden transition 
from a floating to a pinned structure. In this respect, 
this transition represents a dynamical analogue of the 
Aubry transition taking place in the static model at Ac. 
The important difference is that, for the dynamical case, 
this transition occurs as a function of time at all values 
A < A.. 
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FIG. 6. N = 144, T = 233/144, A = 0.05, Po = 0.20. 
(a) Time evolution of log(p^ — Co), (b) CM momentum P. 
Note that the CM momentum has begun its decay before the 
deviation of from its constant value occurs. 

Before concluding this section it is important to discuss 
the relation between this "dynamical Aubry" transition 
described above and the onset of friction. The onset of 
friction is driven by the coupling of the CM to the mode 
with the modulation wavevector q or its harmonics and 
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consists in a special kind of parametric resonances in- 
volving more than one phonon and where the time de- 
pendent driving terms are themselves in resonance 
The appearance of Umklapp requires instead a phonon 
with finite amplitude and very special (zone boundary) 
wavevector. We could say that the last process is more 
difficult to achieve. In Fig. ||, it can be seen that the 
GAM stays conserved, even after the CM momentum has 
begun its decay. This means that there is an interval of 
time in which the first mechanism is active whereas the 
second has not yet taken place. We can identify two 
times: one which characterizes the onset of friction and 
one which describes the pinning of the system. In this 
respect the dynamical model is much richer than the stat- 
ical one. At the Aubry transition, the appearance of a 
static friction occurs by definition at the same value of A 
at which the system gets pinned. 



analytical results concerning the existence of a new in- 
tegral of motion which represents a generalized angular 
momentum related to a phase invariance in incommen- 
surate systems, and we have confirmed this finding by 
means of numerical simulations. We have also shown 
that, during the dynamics, a breakdown of the conserva- 
tion of the GAM occurs at a well-defined time, signalling 
a dynamical transition from a floating phase to a locked 
one. We have been able to prove that this transition is 
related to the apperance of Umklapp processes, caused 
by non linear couplings of the system. We are currently 
trying to further characterize the nature (order) of the 
transition of the dynamical model. We have furthermore 
shown that the onset of friction and the pinning of the 
system are related but distinct phenomena occurring in 
general at different times, which we have been able to 
identify. 



IV. CONCLUSIONS 
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APPENDIX: 



We are going to provide in this appendix an explicit derivation that =^ 0. Let us consider the first term on the 
right hand side of Eq. (13): 



nQX—fiqXnq 
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In order to simplify the notation, we will adopt from now on the following convention: nq = k. Take the derivative 
ofEq. (|Al|): 



(A2) 



We can immediately see that for symmetry reasons the second term cancels. Let us now take the equation for Xk as 
follows from the Euler-Lagrange equations: 
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and the equation for the CM motion: 
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Let us insert Eq. (A3) in Eq. ([A2|). We get: 
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(Note that the first term in Eq. (A3) cancels for the same symmetry reasons given above). We have for simphcity 
exphcitly written down only the first part of the expression in square brackets, since we treat the second part exactly 
in the same way. Rearranging the delta function and applying the following symmetry transformation: 



Eq. (A5) becomes: 
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Because there is no preferential order in the K-summation, the following equality holds: 
There are (m + 1) possibilities, thus we can make the substitution: 

K + Kl + • • • + K„i 

K = 

TO + 1 



Therefore, Eq. (A5) becomes 
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Now, under the assumption that there is no Umklapp ki + • • • = q and can be taken outside the summation. 
Hence, we get: 
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The expression in parenthesis is precisely Eq. (A4) for Q niultiphed by q/27T. Hence, we find: 
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